Call to adopt a nominal set of astrophysical parameters and constants to 
improve the accuracy of fundamental physical properties of stars. 
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ABSTRACT 

The increasing precision of astronomical observations of stars and stellar systems is gradually getting 
to a level where the use of slightly different values of the solar mass, radius and luminosity, as well 
as different values of fundamental physical constants, can lead to measurable systematic differences in 
the determination of basic physical properties. An equivalent issue with an inconsistent value of the 
speed of light was resolved by adopting a nominal value that is constant and has no error associated with 
it. Analogously, we suggest that the systematic error in stellar parameters may be eliminated by: (1) 
replacing the solar radius R e and luminosity L Q by the nominal values that are by definition exact and 
expressed in SI units: = 6.95508 x 10 s m, and l£g = 3.846 x 10 26 W; (2) computing stellar masses 
in terms of M by noting that the measurement error of the product GM Q is 5 orders of magnitude smaller 
than the error in G; (3) computing stellar masses and temperatures in SI units by using the derived values 



At, 



2010 



= 1.988547 x 10 Ju kg and = 5779.57 K; and (4) clearly stating the reference for the values 



of the fundamental physical constants used. We discuss the need and demonstrate the advantages of such 
a paradigm shift. 



Subject headings: Stars: fundamental parameters 
1. INTRODUCTION AND MOTIVATION 

Notable improvements in the precision of as- 
tronomical observations start to challenge our cur- 
rent models and understanding of the physical pro- 
cesses in stars and stellar systems. Until recently, 
the accuracy of fundamental stellar parameters de- 
rived from observations was limited by the stochas- 
tic uncertainty of observations, but w ith new space- 



borne instrumen ts such as MOST (Rucihsk i et al. 



2003b. CoRoT dAuverene et al.1 |2009|) and Kepler 



dBorucki et al.1120101) . systematic effects due to model 
inadequacies are starting to play an important role. 
One of the increasingly important deficiencies is 



the use of the solar mass M , solar radius R Q , and 
solar luminosity L as units in which the fundamen- 
tal stellar properties, deduced via fitting the observa- 
tions by appropriate models, are usually expressed. 
Such stellar properties inevitably depend on the par- 
ticular values of the fundamental solar characteris- 
tics adopted by different researchers in their studies. 
This inconsistency is rooted deeply in the literature 
because the effect of discrepant parameters was of- 
ten well within the systematic and stochastic uncer- 
tainties of observations. However, in certain fields, 
such as binary star research, the associated artifacts 
become increasingly important, especially with the ac- 
cumulation of precise observations of systems with 
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Table 1 : The values of 1 solar mass in the last 35 years. 



M [lO 3O kg] 


Source: 


1.9891 


1 (1976 1 AU constants) 


1.991 +0.002 


2 (Handbook of Chem. Phys.) 


1.98892 ±0.00025 


3 (Astron. Almanac 1994) 


1.9884 ±0.0002 


1 (Astron. Almanac 2011) 



Column "Source" references: 

1... http://asa.usno.nnvy.mil/SecL/Comstants.html 2... Weast & Astie 
i 1980); 3... Astronomical Almanac for the year 1994, U.S. Government 
Printing Office, Washington and Her Majesty's Stationery Office, London 
(1993) 



longer orbital periods. They actually remain a funda- 
mental issue across all fields because eclipsing bina- 
ries (EBs) are the most accurate and commonly used 
calibrators for the masses and radii of single stars 
dTorres et al. 2010h . A typical accuracy of stellar pa- 
rameters for well-studied EBs is claimed to be bet- 
ter than ~2% (And ersen! 1991 ). The modern methods 
for spectroscopic analysis, such as cross-correlation 



dSimkinl [19741 iDa Costa et all 1 19771: iTonrv & Davis 
1979b, broadening functions JRucinskil 119921119981 
and d isentangling ( Simon & Sturml Il994t lHadrav; 



19951 119971: llliiic et all 120041). as well as improved 
light-curve synthesis codes dWilson & Deyinney|l97ll: 
IWilsonll 1 979l 1 1 994 12OO8I: IPrsa & Zwitterll2005N2006l) 
can push these limits even further. The systematic er- 
ror due to inconsistently used solar parameters hence 
propagates to the derived M-R-L-T calibrations, and 
has a notably adverse effect on the computed absolute 
scales and distances. The situation, however, can be 
significantly improved if we are willing to part with 
the current values of canonical solar parameters and 
replace them with the nominal values. 

This issue is not novel in astronomy and a prece- 
dent has been set with the exact value of the speed of 
light. In 1975, Resolution 2 of the 15th Conference 
Generale des Poids et Mesures (CGPM) proposed the 
value of the speed of light in vacuum to be c - 
299 792 458 ms~'. The value was chosen considering 
the excellent agreement to 5c I c ~ 4 x 10~ 9 among dif- 
ferent measurement methods across all wavelengths. 
In 1983, an argument was made that the unchanging 
speed of light needs to be maintained, notably for as- 



tronomy and geodesy, and this was formalized by Res- 
olution 1 of the 17th CGPM that defines the speed of 
light to be the value recommended in 1975. Conse- 
quently, 1 m was redefined to be the distance traveled 
by light in vacuum in 1/299792458 s. This way the 
nominal value of 1 m depends on a fundamental natu- 
ral constant rather than the other way around, as was 
historically the cas^H 

The re-definition of the bolometric magnitude fol- 
lowed at the 23 rd General Assembly of the Interna- 
tional Astronomical Union, held in Kyoto in 1997 Au- 
gust, with adopting a resolution that specifies that 
the zero point of the bolometric magnitude scale will 
no longer be defined through the bolometric luminos- 
ity of the Sun, but rather by defining that M DO \ = 
0.0 mag corresponds to the bolometric luminosity 
Lboi = 3.055 x 10 28 W. This introduces an absolute 
scale of bolometric magnitudes, M\, Q \ = 71.2125 mag- 
2.5 log L, where the bolometric luminosity L is given 
in watts. The convenience of comparing L to the solar 
value remains, since a new definition complies with the 
most often quoted value of the solar bolometric mag- 
nitude: L Q = 3.846 x 10 26 W => M bol>0 = 4.75 mag. 

Another important deficiency is the use of inconsis- 
tent (often outdated) values of the fundamental phys- 
ical constants without providing a reference to the 
source of the value used. The derived constants cannot 
be made exact since they observationally depend on 
fundamental SI units. The universal constant of gravi- 
tation G, for example, is defined as the proportionality 
constant in Newton's law and is one of the mo st dif- 
ficult constants to measure to a high accuracy (Gillies 



1997). The currently recommended 2010 value of the 
constant by the Committee on Data for Science and 
Technology (CODATAflis G = (6.67384 ± 0.00080) x 
10~ n m 3 kg _I s -2 . Since it is derived, this value is sub- 
ject to change in the future as more precise measure- 
ments become available. The solution to this problem 
is to commit to the use of constants set forth by the 
IAU or COD ATA and meticulously provide a reference 
to the used value. 

In this article we quantify the systematic effect of 
inconsistent values of various parameters and propose 
to adopt the nominal values for the mass, radius and 
luminosity of the Sun. We further propose to unify 
the astrophysical constants across modeling tools and 



'See http://www.bipm.org/en/CGPM/db/ 17/ 1 
2 suggested by Dr. Roger Cayrel from Paris. 
3 http://physics. nist.gov/cuu/Constants/index. html 
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Table 2: The value of mean solar radius R G adopted 
from various sources and the corresponding solar ef- 
fective temperature. The (almost negligible) errors in 
values of the effective temperature are given solely by 
the error in the determination of the Stefan-Boltzmann 
constant <x. 



Re 


T e ff 


Source 


(km) 


(K) 




696200.0 


5776.702 + 0.005 


A 


695990.0 


5777.573 + 0.005 


B 


695508.0 


5779.575 + 0.005 


C 


695835.3 


5778.215 + 0.005 


D 


695833.1 


5778.224 + 0.005 


E 


695770.0 


5778.486 + 0.005 


F 


695658.0 


5778.952 + 0.005 


G 



Abbreviations in column "Source": 

A... IAU system of constants 1976; B... Allen 3rd edition: lAllenl 119761) ; C... 
Allen 4th edition: Cox ( 2000), adopted from Brown & Christensen-Dalseaard 
i 1998); D... this paper: a sinusoidal fit to all Cote d'Azur data 1975-1998; 
E... this paper: a sinusoidal fit to good Cote d'Azur data 1975-1998; F... 



Tripathy & Arrti3 fT999l) ; G... IHaberreiter et al. 12008). 



make a strong effort to keep the values of the de- 
rived constants up to date with the IAU, CGPM and 
CODATA resolutions. An excellent example is a re- 
cent review on ac curate stellar masses and radii by 
Torres et al. who quote values of all the rele- 

vant constants used. 



2. QUANTIFICATION OF THE EFFECT 

The r ecently publishe d Kepler Eclipsing B inary 
Catalog jPrsa et al] l201ll; Isiawson et ail l201ll) con- 
tains over 2200 EBs in the 105 deg 2 field of view. 
The period distribution due to Kepler's uninterrupted 
baseline does not suffer from any significant selec- 
tion effects toward longer periods (~10 2 days) and a 
significant number of sources with P orD > 50 d are 
emerging from the sample. Let us consider as an ex- 
ample a binary system with two 1 M components 
in a circular orbit with the orbital period of 200 d. 
The separation between the components can be read- 
ily computed from Kepler's third law. If we adopt the 
above quoted value of G and compute the separation 
based on the 1 M© values listed in Table Q] we will 
arrive at the following values (in 10 9 m): 126.1624 + 



0.0050, 126.2026 + 0.0050, 126.1586 + 0.0050, and 
126. 1476 + 0.0050. This accounts for the relative error 
of 4.0 x 1 -5 , or an absolute error of 5 .0 x 1 6 m. A sin- 
gle event in Kepler long-cadence data (30-minute ex- 
posure) can be timed to ~6 minutes; a P 01D = 200-day 
binary will have 1 8 events (nine primary and nine sec- 
ondary eclipses) over the 5 yr mission lifetime, which 
will reduce the timing error to ~1.4 minutes, or in rel- 
ative terms to 5 x 10~ 6 . This is an order of magnitude 
smaller than the effect of using inconsistent solar-mass 
values, and will be the cause of significant systematics. 
If short-cadence data (1 -minute exposures) are avail- 
able for the given target, the timing is improved by 
another order of magnitude, and the systematics will 
overpower the stochastic error. 

The same is true for single stars as well. Consider 
a giant star with a radius of 30 R and an equatorial 
rotational velocity of 5 km s _1 . The period of rota- 
tion of such a star would then be 303.801 days for the 
IAU 1976 va lue of the solar radius but 303.473 days 
for the recent iBrown & Christensen-Dalsga ard (1998) 
value (see Table 2). Such a difference is readily de- 
tectable after a few rotational periods covered by rele- 
vant observations. 

Also worth considering is the dependence of the 
effective temperature on the value of stellar radius 
through L = AjiR\crT^ s , where <x is the Stefan- 
Boltzmann constant. Table [2] lists several values of 
the solar radius, either recommended in various com- 
pilations or derived from recent measurements. For 
each value we give the corresponding effective tem- 
perature calculated for the solar luminosity value of 
L = 3.846 x 10 26 W and cr = 5.670373(21) x 
10~ 8 WirT 2 K~ 4 . The systematic differences in T e ff 
values are larger than the propagated errors. 

3. THE SOLAR UNITS 

It is customary in stellar research to express the ba- 
sic physical properties of the stars, such as their lu- 
minosity, mass or radius, in solar units. While conve- 
nient, this is somewhat unfortunate, for the following 
reasons: 

1. the values of the solar luminosity, mass and 
radius are subject to continuous improvement 
thanks to the increasingly precise observational 
techniques; and 

2. solar luminosity and radius vary measurably 
with the solar cycle and perhaps also on other, 
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Table 3 

Select examples of the commonly used equations utilizing the nominal solar values proposed by this paper, and 

constants recommended by the iau. 



Select example: 



Equation: 



Comments: 



Kepler's third law: 

Surface gravity: 

Stellar bolometric magnitude: 
Parallactic radius: 

Angular orbit diameter: 
Equatorial rotational velocity: 

Mass functions: 



A[m] = 2927699260.629(74)/' 7 (M, + M 2 )3 
A[AU] = 0.01957046077547(49)/^ (M, + M 2 )3 

A[Rg] = 4.20944009361(1 1)/>§(M, + M 2 )3 
logg [cgs] = 4.438307381330(33) + logM- 21ogR 

M b0 , = 42.3689588(40) -5 log/?- 101ogr dr 
R = 107.5457584245(22)8 p" 1 

a = (215.0915168490(43))-'/* p 
V = 50.57877 RP^. 



Mi sin 3 i 
M 2 sin 3 a 
MM) 
MM) 



1.036149050206(78) x IPT 1 K 2 (K t + K 2 ) 2 P(l - e 2 ) 1 ' 5 
1.036149050206(78) x 10~ 7 K,(K, + K 2 ) 2 P{\ - e 2 ) 1 5 
1.036149050206(78) x IQ' 1 K]P(\ - c 2 ) 1 5 
1.036149050206(78) x 10" 7 AT| J P(1 - c 2 ) 1 5 



A: semi-major axis in m, AU, and M\ and M 2 : point masses in M 2010 ; P: period in mean 
solar days. 



M: stellar mass in M™°\ R: stellar radius in Kg. With the use of the nominal values of the 
solar mass and radius, the error only depends on the error of the GM product constant and 
the conversion formula becomes very accurate, highly exceeeding the current precision of our 
knowledge of stellar masses and radii. 

R: stellar radius in ; T c n: effective temperature in K. The error in the conversion constant 
is due to the uncertainty in the Stefan-Boltzmann constant cr. 

R: radius in R^J; 6: angular stellar diameter in arcsec; p: stellar parallax in arcsec. 
The constant is computed from the recommended value of 1 AU and 1 pc by IAU 2009 
(149, 597, 870, 700(3) m and 3.085677581503(62) X 10 16 m, respectively). 

a: angular orbit diameter in arcsec; A: semi-major axis in au; p: stellar parallax in arcsec. The 
same constants were used as above. 

V: equatorial rotational velocity in km s~' ; R: stellar radius in Kg ; P M : rotational period in 
mean solar days. Note that the constant in the above equation has no uncertainty thanks to our 
use on the nominal value of the solar radius. 



Mi : stellar mass in M 010 ; i: inclination in deg; K\ 2 : RV semi-amplitudes in km s~' ; P: 
orbital period in mean solar days; e: eccentricity. The uncertainty stems from the uncertainty 
in the GM Q product. 



Projected orbital sizes: 



a\ sin( 
a 2 sin (' 
A sin i 



0.019771 \\2K 2 P(\ -e 2 ) 05 
0.019771142(Kj + K 2 )/>(!- e 2 f 



Solar effective temperature: 



T c«. = 5779.5747(54) 



a\X- distance of the orbiting component from the center of mass, in H^; A = a\ + a 2 : semi- 
major axis in Jt^; i: inclination in degrees; K\ 2 : RV semi-amplitudes in kms _I ; P: orbital 
period in days; e: orbital eccentricity. Note that there is no error in the transformation constant. 

T e j is the effective temparature of the Sun for the nominal value of the solar radius 



even shorter time scales (see, e.g. JSelhorst et al. 
2004). Due to effects such as mass loss via stel 



lar wind or mass gain due to infalling material, 
the change of the solar mass may become mea- 
surable on longer timescales. 

At the same time, some important quantities, such 
as the projected velocity of the stellar rotation, are 
measured in absolute units (km s _1 in this particular 
case) that depend directly on the adopted values of the 
fundamental solar properties. 

We suggest a simple solution to the problem of in- 
consistent use of various different values of the solar 
units - one that follows the precedent for the speed of 
light and bolometric magnitude. Let us deprecate the 
use of the current values of R Q and L with all their un- 
certainties and time variations, and replace them with 
the exact, nominal units 7?q and Xq, where their nu- 
merical values are chosen to be close to the recently 
derived (or adopted) values: the sola r radius from 
Brown & Christensen-Dalsgaardl (119981) . also recom- 
mended in the 4 th edition of Allen's Astrophysical 



Quantities dCox 2000), and the solar luminosity to 
comply with the above-mentioned IAU resolution on 
bolometric magnitude. 

Once the values of R Q and L are made nominal, the 
relationship between the nominal values 9?§ and X@, 
and the effective temperature 7~g 010 can be determined 
readily from: 

Xg = 4 ^V(r| 010 ) 4 ,where 
a = 5.670373(21) x 10~ 8 WnT 2 KT 4 . 

Note that the small error in the determination of the 
Stefan-Boltzmann constant <x propagates to negligibly 
small errors in the determination of the effective tem- 
perature. 

Unfortunately, the same cannot be proposed for 
the solar mass: while the universal gravitational con- 
stant G = (6.67384 ± 0.00080) x 10" 11 m 3 kg~V 2 
is one of the least precisely determined funda- 
mental constants in nature, the product GM G = 
1.32712442099( 10) x 10 20 m 3 s" 2 is determined much 



more precisely ( Petit et al. 2010b . All physical pa 



rameters that can be expressed in terms of GM Q can 
thus be determined much more accurately than the 
ones depending on M alone. Table [3] lists several 
imporant ones, most notably stellar masses in terms 
of solar mass, and the absolute scale of the system 
in SI units. In particular, the mass of a binary star 



component j ( j - 1 , 2) can be expressed as GMj = 
[K 3 -jP(Ki + K 2 ) 2 ]l{2n sin 3 i), where Ki and K 2 are ra- 
dial velocity semi-amplitudes in km s _1 , and i is the in- 
clination. If we divide this expression by GM G , we ob- 
tain Mj/M Q in terms of GM G , which is much more pre- 
cise than computing it with the more uncertain value of 
G. However, both solar and stellar masses in SI units 
will still be limited by the measurement error in G. If 
solar mass were to be made nominal, the precision in 
Mjl would be degraded by 5 orders of magni- 
tude. To make the distinction clear, we use the des- 
ignation N to represent nominal parameters, and 2010 
to denote parameters that are derived from the latest 
values of GM Q and <x. 

The proposed nominal and derived parameters to be 
used in computations are: 



1 = 6.95508 x 10 s m 



\3 
1 Xq 

1 GM 2010 

y2010 

-2010 



I Ml 



1 T, 



© 



3.846 x 10 26 W 

1.32712442099(10) x 10 20 mV 2 
1.988547 x 10 30 kg 
5779.57 K 



4. SELECT EXAMPLES 

There are a number of frequently used formulae 
where numerical constants are affected by the adopted 
values of the solar mass, radius and luminosity. We 
list some of them for the suggested nominal values 
in Table [3] In all examples, the propagated errors of 
the fundamental physical constants used are given in 
parentheses and denote the uncertainty of the last two 
significant digits. The expressions for the semimajor 
axis A from the third Kepler law are given for the value 
of A in km, AU, and <Rg. 

5. CONCLUSIONS 

We presented the case for obsoleting the use of R Q , 
L Q and M Q as units and replacing two of them with 
their nominal counterparts ^ and Xq that are exact. 
In all determinations of stellar masses expressed in 
M Q , the product GM^ m , known to a very high pre- 
cision, should be used. For the conversion of the thus 
derived stellar masses to SI units, Al 2 , 010 can be used. 
Once the precision of the universal gravitational con- 
stant G is improved, it will be beneficial to deprecate 
the actual value of the solar mass and replace it by an 
exact, nominal value Mq. 
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All the recommendations presented in this article 
would reduce the systematics that stem from using dis- 
crepant values of fundamental properties of the Sun. 
This is a consequence of the true variations of these 
values due to intrinsic effects such as the magnetic so- 
lar cycle, as well as of a steady improvement in their 
determination via more precise observations. We fur- 
ther implore the community to use fundamental phys- 
ical constants recommended by the IAU and/or CO- 
DATA and call for meticulous referencing of the used 
sources. 
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SC-1041. We acknowledge the use of the electronic 
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